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Abstract. In this paper, we investigate the epidemic spreading for SIR model in 
weighted scale-free networks with nonlinear infectivity, where the transmission rate in 
our analytical model is weighted. Concretely, we introduce the infectivity exponent 
a and the weight exponent /3 into the analytical SIR model, then examine the 
combination effects of a and (3 on the epidemic threshold and phase transition. We 
show that one can adjust the values of a and f3 to rebuild the epidemic threshold to 
a finite value, and it is observed that the steady epidemic prevalence R grows in an 
exponential form in the early stage, then follows hierarchical dynamics. Furthermore, 
we find a is more sensitive than /3 in the transformation of the epidemic threshold and 
epidemic prevalence, which might deliver some useful information or new insights in 
the epidemic spreading and the correlative immunization schemes. 
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1. Introduction 

There has been a long history for the research of epidemic spreading [1]. And in 
the general case, the epidemic system can be represented as a network where nodes 
stand for individuals and an edge connecting two nodes denotes the interaction between 
individuals. In the past, researchers mainly focused the disease transmission study on 
the conventional networks [2, 3] such as lattices, regular tree, and ER random graph. 
Since late 1990s, scientists have presented a series of statistical complex topological 
characteristics [4, 5, 6, 7] such as the small-world (SM) phenomenon [8] and scale-free 
(SF) property [9] by investigating many real networks including the internet [10], the 
www [11], the scientific web [13], the protein networks [14] and so on. Subsequently, the 
studies of dynamical processes on complex networks also have attracted lots of interests 
with various subjects [15, 16], and as one of the typical dynamical processes built on 
complex networks, epidemic spreading has been investigating intensively once more. 

The basic conceptual tools in understanding the epidemic spreading and the related 
effective strategies for epidemic controlling should be the epidemiological models [17, 18]. 
Among the numerous possible models, the most investigated and classical models are SI 
model [19, 20], SIS model [21, 22, 23, 24, 25] and SIR model [26, 27, 28, 29], which can 
approximately describe the spreading of real viruses such as HIV, encephalitis, influenza 
virus in biological networks; computer virus, trash mail in technological networks; and 
even gossip in social networks. The most valuable result in standard SIR model (or SIS 
model) is that the critical threshold (of transmission rate) vanishes for the scale-free 
networks in the limit of infinite network size [28]. 

In consideration of epidemic spreading in real cases, there yet has been some 
inappropriate assumption in the details of the standard SIR model. We know that 
in the classical SIR model the transmission rate A is a constant, but in the real world, 
A should be different among individuals. Based on this assumption, in reference [30], 
Jaewook Joo et al. proposed the effective transmission rate that they introduced a 
effective coefficient C{k,l) based on standard transmission rate A for the edge {k,l); 
similarly, in reference [31], Ronen Olinky et al. also studied the effectiveness of the 
transmission rate, and in their work, the transmission rate is XA{k) where A{k) means 
the probability that a susceptible node actually acquires the epidemic through an edge 
connected an infected node with degree k. In these previous studies, the transmission 
rate on a given edge is just treated as a function with degrees of the two connecting 
nodes, which will induce the transmission rates of two opposite directions on the same 
edge are symmetrical. And, their analytical methods and results which are based on 
the stationary state of the SIS model may be not valid in the SIR model. 

Thus, in order to make the transmission rate accord with the realistic cases much 
more, we take into account the effects of the weights of edges and the strengths of nodes 
which are of great importance measures in the weighted networks [33, 34, 35]. And 
indeed, the weight (or the strength) is one of the most important indications in lots 
of real networks, for example, in social networks it can represent the intimacy between 
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individuals; in the internet the weight can imply the knowledge of its traffic flow or 
the bandwidths of routers [7]; in the world-wide airport networks it can evaluate the 
importance of a airport [5, 33], and so on [34, 35]. Particularly, for epidemic spreading, 
the weight can indicate the extent of frequency of the contacting of two nodes in scale- 
free networks, the larger the weight is, the more intensively the two nodes communicate, 
at the same time, the more possible a susceptible individual will be infected through 
the edge where the transmission rate is larger. 

On the other hand, in the classical SIR model, each infected individual can establish 
contacts with all his/her acquaintances (neighbors) within one time step, that is to say, 
each infected node's infectivity equals its degree. But in the real case, a individual can't 
contact all his intimate friends, particularly when he is a patient. In reference [38], 
the infectivity is assumed as a constant A, which means each infected individual will 
generate A contacts at each time step. Recently, Fu et al. proposed a piecewise linear 
infectivity [39], which means: if the degree /c of a node is small, its infectivity is ak; 
otherwise its infectivity is A as a saturated value when k is beyond a constant A/ a. 
Both the constant or the piecewise linear method, the heterogeneous infectivity of nodes 
with different degrees is not considered as adequately as possible in scale-free networks, 
that is to say there may be some nodes with different degrees which have the same 
infectivity, and there will be a large number of such nodes if the constant A is assigned 
irrelevantly or the size of underlying networks is infinite. So, in order to solve these 
problems, we introduce the nonlinear infectivity, namely, a infectivity exponent a will 
be introduced to take control of the number of contacts that a infected node generates 
within one time step, and the a is between and 1 which is convenient to adjust for 
different scale-free networks. 

In this paper, we present the modified SIR model where the infectivity exponent 
a and the weight exponent (3 are added; based on the modified model, the dynamical 
differential equations for epidemic spreading is proposed. We parse the equations to 
investigate the threshold behavior and propagation behavior for epidemic spreading; 
and the analytical results we obtain is verified by the necessary numerical stimulations. 
We show that one can adjust the values of a and (3 to rebuild the epidemic threshold to 
a nonzero finite value for different networks, which can prohibit or delay the epidemic 
outbreaks to some extent. And we find a is more sensitive than f] in the transformation 
of the epidemic threshold and epidemic prevalence, which indicates the intrinsic factor 
(the infectivity exponent a) take more responsibility than the extrinsic factor (the weight 
exponent j3) for the epidemic outbreaks in large scale-free networks. 

2. Standard SIR model 

Epidemic modeling has a history of researching, and mathematicians also put forward 
many epidemic models [2, 3, 18]. In the domain of complex networks, SIR model [26, 27, 
28, 29] is one of the most investigated and classical epidemic models. In the standard 
SIR model, individuals can be divided into three classes depending on their states: 
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susceptible (healthy), infected and removed (immunized or dead). In order to take 
into account the heterogeneity induced by the presence of nodes with different degrees, 
we use Sk{t), Ik{t), Rkit) to denote the densities of susceptible, infected and removed 
individuals with degree k at time t, respectively. And these variables are connected 
by means of the normalization: Skit) + Ikif) + Rkit) = 1. The global quantities such 
as the (average) epidemic prevalence are therefore expressed by an average over the 
various degree classes, i.e., -R(t) = J2k^i^)^k{t)- For the standard SIR model, the 
epidemic evolves by the following rules: at each time step, a susceptible individual 
acquires the infection at the transmission rate A in one contact with any neighboring 
infected individual, which means if a susceptible individual has a edge connecting a 
infected individual, the disease will transmitted to the susceptible one through the edge 
with a specific probability A. On the other hand, the infected ones will recover and 
become immune (can't be infected any more) with rate fi, one can set /i = 1 without 
loss of generality. 

For a comparison, firstly we review some classical results from reference [28], where 
Moreno et al. used the mean field theory to describe the dynamical differential equations 
of SIR model as follows: 

Rk{t))J2Pik'/k)Ik^{t), (1) 
k' 

hit)- Rkit)) J2Pik'/k)Ik'{t), (2) 
k' 

(3) 

where P(fc'/fc) denotes the conditional probability for a node with degree k to connect 
a node with degree k' . In the uncorrelated case, Moreno et al. obtained the epidemic 
threshold: Ac = {k)/{k'^), which implies the absence of the epidemic threshold in a 
wide range of scale-free networks ((A;^) — > cx), Ac — >■ 0). This result is a bad message 
for epidemic controlling and preventing, since the epidemic will prevail in many real 
networks with any nonzero value of transmission rate A. 



dSkit) 
dt 

dikit) 
dt 

dRkit) 
dt 



-Xkil-Ikit) 
-Ikit) + Xkil- 
- hit), 



3. SIR model in weighted network 

In this section, we will give a detailed investigation about the modified SIR model into 
which we introduce the weighted transmission rate and nonlinear infectivity. The results 
we obtain might deliver some useful information for the epidemiology. And for a better 
analysis, we firstly describe the general differential equations for SIR model based on 
the mean field theory, as follows: 

^ = -kil- hit)- Rkit)) J2Pik'/k)h'it)^Xk'k, (4) 

k' 
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^ = + HI - hit) - Rk{t))J2P(k'/k)h'it)^X,',, (5) 

k' 

dRk{t) 

= /^^H^)' (6) 

where Sk, Ik, Rk have the same meaning with the standard SIR model (see section 2); 
and ^{k'), denote the infectivity of nodes with degree fc' and the transmission rate 
from nodes with degree fc' to nodes with degree fc, respectively. 

3.1. The model 

Different from the previous studies, in this paper, we mainly focus the SIR model on the 
weighted networks. Among varieties of weighted patterns in complex networks, making 
use of nodes's degrees to express the weights of edges is very important, namely, the 
weight between two nodes with degree k and k' may represent as a function of their 
degrees [33, 34, 35], i.e., w^^' = wo{kk')'^, where the basic parameter wq and the exponent 
(3 depend on the particular complex networks (e.g., in the E.coli matabolic network (3 
= 0.5; in the US airport network (USAN) f3 = 0.8 [36]; in the scientist collaboration 
networks (SCN) /? = [33]). Noteworthily, the weight Wf^i^' belongs to an edge, similarly, 
a node (with degree k) also can be measured by weights, i.e., the strength of a node (with 
degree k), which can be obtained by summing the weights of the links that connected 
to it, i.e., Nk = kTij^i P{kl /k)wj^^i , where Nk is the strength of a node with degree k. In 
this paper, for simplicity, we focus on uncorrelated (also called non- assort at ive mixing) 
networks where the conditional probability satisfies P{k' /k) = k' P(k')/ (k) [37]. Thus, 
one can obtain A^^ = WQ{k^'^^)k^^^ / (k). 

Here, for each node with degree k we fixed a total transmission rate which is given 
by A/c, and the transmission rate on the edge from the /c-degree node to /c'-degree node, 
will be redistributed by the proportion of the k-degree node's strength that the edge's 
weight accounts for, that's to say the A^^' can be defined as follows: 

V = AA:^, (7) 

from which we know the more proportion of Nk that the weight Wj./^' of an edge accounts 
for, the more possible the disease will transmit through the edge. In the uncorrelated 
case, one can obtain A^^' = \k'^{k) / {k^^^^). Moreover, the reasonable total probability 
that a susceptible node with degree k will be infected at time step t is given by 
-'- ~ Y\\tk' (^NiDS{t)0- ~ ^kk')y where NIDS(t) denotes the degree sequence of neighboring 
infected nodes that connect to the susceptible node with degree k at time step t. 

On the other hand, from the general differential equations of SIR model (equations 
(4) - (6)), we know (f{k) denotes the infectivity of nodes with degree k, and here, in the 
present model we define it as follows: 

ip{k) = k'',0<a<l, (8) 



Epidemic spreading with nonlinear infectivity in weighted scale-free networks 



6 



which is to say, each infected individual can estabhsh contacts with its fc" neighbors 
within one time step. The exponent a will dominate the infectivity among nodes with 
different degrees. Since < a < 1, it can be adjusted to make the contacts fall on a 
more realistic range. And the node's infectivity will grow nonlinearly with the increasing 
degree k. 

We take the simplified expressions of (p{k) and A^.^y into the equations (4) - (6) with 
= 1 (without lack of generality), we obtain as follows: 

dSkit) Xk^+P 



dt 



Sk{mt), (9) 



^ = -h{t) + ^s,mt), (10) 
^ = m (11) 

where 6{t) = Y2k k°'P{k)Ik{t)- The above equations combined with the initial conditions 
R^{0) = 0, 4(0) = and 5^(0) = 1 - 4(0) - Rk{0) = 1 - /°. And in the general case 
is very small, then we can obtain Sk{0) ~ 1. Under this approximation, equation (9) 
can be directly integrated, as follows: 

Skit) = e~(^^^*\ (12) 

where 0(t) = 9{t)dt = ^^j. k"'P{k)Rk{t), and in the last equality we have made use of 
equation (11). 

3.2. Threshold behavior 

In order to obtain some material results for the epidemic threshold and the average 
epidemic prevalence, firstly we compute the time derivative of the magnitude [28]: 

^=J2k-p(k)[l-R,{t)-S,{t)] 

k 

= (O-0(t)-^A;"P(A;)5fc 



= (r) - (Pit) - J2 k''P{k)e (13) 

k 

For the general P{k) distribution, equation (13) can not be solved in a closed form. 
However, we can still obtain some useful results in the steady state of the epidemics. 
Since in the steady stage with sufficiently large t, we have that 4(oo) = 4 = and 
consequently limf^oo d0(t)/dt=O, then one can get the self-consistent equation for (f) 
from equation (13) as follows: 

= (A;") - ^ A;"P(A;)e (14) 
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The value 0=0 is always a (trivial) solution. Then we compute the second order 
derivative of the rhs of equation (14) for </>, and note that 

— ((A;") - J2 k'^P{k)e~T^)^) < 0, (15) 

we can see the rhs of equation (14) is a convex function, therefore, a nontrivial solution 
of equation (14) exists only if the condition 

Hk'^) .J2k''me~^\=o > 1, (16) 

can be satisfied. This relation implies 

Y.m^^ = x^>i, (17) 



k 



and the above inequation defines the epidemic threshold: 

A - ^^'^'^ (18) 

Below which the average epidemic prevalence {R{t)) will finally be approximatively 
null, and above which it will attain a finite value. One can see if a = l,/5 = 0, then 
Ac = {k)/{k'^), which induces the absence of the epidemic threshold in a wide range of 
scale-free networks [32]. And if a + f3 = 0, the threshold will be a finite value given by 
Ac = {k^~^^) / (k) = {k^~"') / (k) > 1/ {k); similarly, if a-\- P = —1, one can also get a finite 
threshold which is Ac = = (A;"") > (l/k). 

Furthermore, we consider the epidemic threshold in the case of general scale-free 
networks of which the degree distribution is P{k) = ck~'^, 2 < 7 < 3, where c is the 
normalization constant. Then, we obtain {k^^'^) = c{{k^'^~^) — {k^^~^)) / {(3 + 2 — 7) 
and (A;°+^+i) = c{{k^,t^-') - (Ct„^+'"''))/(« + /5 + 2 - 7), where fc^ax {k^n) denotes 
the largest (smallest) degree in the underlying networks. Substituting into equation 
(18), one can rehandle the epidemic threshold as follows: 

' /3+2-7 k'^^i^'-'-k':^^'-'' 

From equation (19), one can see that the infinite of the largest degree (/Cmax 00 or 
equally N 00, since /cmax oc A^^/('>'~^) [12]) will make the epidemic threshold Ac tends 
towards zero if 7 < a + (3 + 2; on the other hand, if 7 > a + f3 + 2, the epidemic threshold 
Ac is approximate to be a finite value, given by Ac = A;~° (a -I-/3 + 2 — 7)/(/5-|-2 — 7). 
Thus, the critical border is 7 = a + /5 + 2. Although for most real networks including 
the internet [10], the www [11], the world-wide airport networks [33] and the scientific 
collaborations networks [13], the topology exponent 7 exists between 2 and 3, which is 
incidental to induce the absence of the epidemic threshold, one can adjust the infectivity 
exponent a and the weight exponent (3 to restore a nonzero threshold for a given 
networks (a fixed value of 7) . 

In order to get a intuitionistic relation with Ac, a and P, we have performed 
numerical simulations for the epidemic threshold. Firstly, on the basis of the simulated 
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Figure 1. (a): The epidemic threshold Ac versus a with the exponent /? = in BA 
networks, (b): The epidemic threshold Ac versus f3 with the exponent a = 1.0 in BA 
networks. 




Figure 2. The epidemic threshold Ac in a 3D-graph, which is made up with various 
a and (3 in BA networks. 



stochastic realizations of SIR model on BA networks [9] of which the theoretical scale- 
free exponent is 3 (7 = 3), we fix /5 = and adjust a between and 1 to show 
the transformation of epidemic threshold Ac- In this case the critical value of a is: 
ac = 7 — /? — 2 = 1, below which Ac is a nonzero finite value. As figure 1(a) displays, 
the value of Ac is greater than 0.01 with the relation a <^ 1, vice versa. Secondly, we fix 
a = 1.0 and adjust j3 between -2 and 2 to show the transformation of Ac. In this case 
the critical value of /5 is: /3c = 7 — a — 2 = 0, below which Ac is a nonzero finite value. 
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X 



Figure 3. The steady epidemic prevalence R versus A for SIR model in BA networks 
with N = 10'', (k) = 6, /3 = 0, and a=1.0, 0.9, 0.8, • • •, 0.2, 0.1 (from top to bottom). 



As figure 1(b) displays, at the point of /3c = 0, Ac — 0.03, and Ac will take on a mucli 
faster cliange in /5 < than the one in /3 > 0, and the finiteness of Ac is apparent when 
P < —0.5. From figure 1, one can see the simulations are consistent with the analytic 
results about critical threshold when we consider the effect of finite scale of the substrate 
work we use. Moreover, it is observed that the decreasing trend of Ac with increasing a 
is much quicker than the one with increasing j3, that is to say, a is more sensitive than /3 
in the transformation of epidemic threshold Ac, which means a is the leading factor for 
the transformation of Ac in the present model. Figure 2 displays the epidemic threshold 
in a 3D-graph, which is made up with various a and /3, and the critical condition for 
a nonzero finite threshold is a + (3 < 1. One can see that Ac is small in the blue area 
where the great mass of data meet the condition that a + /5 > 1, which is the condition 
of threshold vanishing. 



3.3. Propagation behavior 

For further investigation of the epidemic dynamics of the present model, we study 
the propagation behavior of the epidemic spreading. Firstly, we investigate the 
average epidemic prevalence in the steady stage of epidemic evolution with different 
combinations of a and (3. From the analysis in section 3, it is easily to conclude that 
/fc = at the epidemic critical point of A, which will induce a quite small Rk, and we 
can approximately get ~ according to the relationship = k°'P{k)Rk. Then, 
expanding the rhs of equation (14) for the small 0, and ignoring the higher-order terms, 
we obtain 



-A 



l(a + 2^+2) 2. 
2 {k^^f^y 



1. 



(20) 



Epidemic spreading with nonlinear infectivity in weighted scale-free networks 



10 




Next, we compute the derivative of equation (20) for A at the critical point, as follows 

2(r+'3+i)3 



As referred above, one can obtain that 



(21) 



R = Y,P{k)Rk = 1 -Y,P{k)e-^'^^'^'^^'^"\ (22) 

k k 

consequently. 
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Figure 5. The steady epidemic prevalence R versus A for SIR model in BA networks 
with N = 10'', (fc) = 6, a + /3 = -1, and from top to bottom: a=0.9, 0.7, 0.5, 0.3, 
0.1; accordingly /3=-1.9, -1.7, -1.5, -1.3, -1.1. The inset shows the susceptible density 
S versus A with the same combination of a and /3. 




Figure 6. The steady epidemic prevalence R versus A for SIR model in BA networks 
with N = lO'', (k) = 6, a + P ^ 0, and from top to bottom: a=:0.8, 0.6, 0.4, 0.2; 
accordingly /3=-0.8, -0.6, -0.4, -0.2. The inset shows the susceptible density S versus 
A with the same combination of a and (3. 



Similarly, in the considering of general scale-free networks as referred above, equation 
(23) can be written as follows: 

2c(a + 2/3 + 3-7) jk^ti^'-' - Ctn"^''")^ ^ 



The obtained results shows, the exponent a will make a primary contribution to the 
velocity of increasing of the steady epidemic prevalence (R) by given the topology of a 
underlying network. Combining the analysis in section 3, for a fixed sum of a and (3, 
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Figure 7. The steady epidemic prevalence R versus A for SIR model in BA networks 
with N = 10^, (fc) = 6, a + /3 = 1, and from top to bottom: a=0.7, 0.6, 0.5, 0.4, 0.3; 
accordingly (3=0.3, 0.4, 0.5, 0.6, 0.7. The inset shows the susceptible density S versus 
A with the same combination of a and (3. 



-I — I — 1 — I — 1 — I — 1 — I — 1 — I — 1 — f 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

A. 



Figure 8. The steady epidemic prevalence R versus A for SIR model in BA networks 
with N = 10'^, (fc) = 6, a + /3 = 2, and from top to bottom: a=0.9, 0.8, 0.7, 0.5, 0.3, 
0.1; accordingly /3=1.1, 1.2, 1.3, 1.5, 1.7, 1.9. The inset shows the susceptible density 
S versus A with the same combination of a and p. 

one can conclude that the more ratio of a to /3 is , the larger the Ac is, and the slowly the 
R grows as A increases (since the most large-scale real networks have the relationship 
a + 1 - 7 < 0). 

For a better understanding of the epidemic propagation behavior, we take numerical 
simulations with various combination of a and (3 on BA networks (7 = 8). Firstly, we 
investigate the impact of a and (3 separately, which are the two particular cases: a = 1.0 
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5 10 15 20 25 30 35 40 45 50 

t 



Figure 9. The temporal epidemic prevalence R{t) versus A for SIR model in BA 
networks with N = 10'*, (k) = 6, a + /? = —1, and from top to bottom: a=0.9, 0.8, 
0.7, • • •, 0.2, 0.1; accordingly (3^-1.9, -1.8, -1.7, • • •, -1.2, -1.1. 

with different (3, and (3 = with different a. Figure 3 displays the effects of a on the 
steady epidemic prevalence R with /3 = 0. As the figure shows, for a same A, one can 
observe the slope of R grows as a increases (see a=0.9, 0.8, ■ ■ ■, 0.2, 0.1), which is 
consistent with the analytical results from equation (24). Figure 4 displays R versus A 
in the case of a = 1.0, and figure 4(a): P=0, 0.5, 1.0, 1.5, 2.5 (from top to bottom); 
figure 4(b): /5=-0.5, -1.0, -1.5, -2.5 (from top to bottom). As shown in figure 4, the 
larger the absolute value of /5 is, the more slowly the R grows. On the other hand, for a 
general case (« 7^ 1 and 7^ 0), according to the critical equation a + /5 + 2 = 7 and the 
algebraic sign of the sum a + (3, we consider four representative combinations of a and 
f3, which are a + (3 = —1, a + (3 = 0, a + (3 = 1 and a + (3 = 2. In each combination, 
there also has been divided into several different configurations by the ratio of a to (3. 
As shown in figures (5 - 8), for a fixed sum of a and (3, one can see the R grows more 
quickly as the ratio increasing, which is consistent with our analytical results that a 
is a more sensitive factor to R. And further more, since a + [3 > 1, one can see the 
epidemic threshold Ac is quite small (tends to zero, considering the effects of finite size) 
in figure 4(a) and figure 8, which is also consistent with the critical condition of Ac; 
otherwise, for a + (3 < 1, Xc becomes to be a nonzero finite value as shown in figure 3, 
figure 4(b), figure 5 and figure 6. For the early stage of A > Ac in figures (3-8), one 
can see the R grows in an exponential form with A increases, then the growth rate will 
take off slowly for A ^ Ac, at last, it will tend to be zero, which means the value of R 
have attained a steady value. 

For investigating the temporal propagation behavior, we simulate the time behavior 
of R{t) for SIR model on BA networks with A = 0.5. As displayed in figure 9 
{a + (3 = -1), figure 10 {a + (3 = 0), figure 11 (a + /3 = 1) and figure 12 {a + (3 = 2), 
one can see that, the prevalence R(t) grows in an exponential form in the early stage, 
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Figure 11. The temporal epidemic prevalence R{t) versus A for SIR model in BA 
networks with N = 10^, (fc) — 6, a + (3 = 1, and from top to bottom: a=0.9, 0.8, 0.7, 
• ■ •, 0.2, 0.1; accordingly l3^0.1, 0.2, 0.3, • ■ •, 0.8, 0.9. 



and then stabilizes in a nonzero value as time goes on. Since A may be smaller than Ac 
when the value of a is much small such as a = 0.1, 0.2, thus if A < Ac, the steady value 
of R{t) (i.e., R) will be quite small, which approximatively equals to the initial density 
of infected nodes; when A > Ac, the prevalence R{t) is higher as the ratio of a to /5 gets 
larger at the same time of t, as the figures display. Moreover, it is observed that the 
steady value of R{t) is smaller in figure 12 compared with the other figures (9 - 11). 
Although the exact solution of R{t) about a and f3 which can demonstrate the difference 
well is difficult to be managed here, from a qualitative perspective, we believe that's 
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Figure 12. The temporal epidemic prevalence R{t) versus A for SIR model in BA 
networks with N — 10'*, (k) — 6, a + P = 2, and from top to bottom: a~0.9, 0.8, 0.7, 
• • •, 0.2, 0.1; accordingly /3=1.1, 1.2, 1.3, • • •, 1.8, 1.9. 



because of the trait of SIR model. In figure 12, a + /3 = 2 which will induce a small 
threshold Ac, thus at the early stage of evolution many susceptible nodes will be infected 
in view of the relation of A Ac, and due to the immune rate we set is unity (/x = 1), 
the old infected nodes will become removed ones that cant't be infected any more at 
the same time. Consequently the optional objects for a infected node will decrease as 
time evolutes. Moreover, there will be many infected nodes surrounded by the removed 
ones. Under this situation, the epidemic spreading will arrive at a equilibrium much 
more quickly, and thus the epidemic prevalence in the steady stage also will be a small 
value as figure 12 displays. 

4. Conclusion and discussion 

To sum up, in this paper, we have investigated the dynamical behavior of SIR model 
with weighted transmission rate and nonlinear infectivity, we present that one can adjust 
the exponent a and /5 to control the epidemic threshold which is absent for the standard 
SIR model in scale-free networks. The critical value just depends on the exponent a and 
f3 for a given topology of networks (a fixed value of 7), and a is more sensitive than f3 
for the transformation of the epidemic threshold and epidemic prevalence, which agrees 
with the numerical simulations very well. And the numerical results of the time behavior 
of R{t) also have been presented, where the remarkable result is that, for a fixed A, the 
smaller threshold Ac will induce a smaller epidemic prevalence at the equilibrium. 

In a way, epidemic spreading can be managed as a reaction-diffusion process [40, 
41, 42], which also has a very close relation with information retrieval, peer-trust and 
influence spreading. The efficient diffusion or inefficient diffusion maybe has its merits in 
various natural and artificial networks. Our work might deliver some useful information 
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or new insights in designing data layout, city layout and network layout for performing 
their best advantages. 
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